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<^ ■ Abstract 

C. 

^ , The effect of quark-gluon vertex dressing on the ground state masses of the u/d-quark pseu- 

^~i ■ 

^^ ' doscalar, vector and axialvector mesons is considered with the Dyson-Schwinger equations. This 

extends the ladder-rainbow Bethe-Salpeter kernel to 2-loops. To render the calculations feasible 
^^ ■ for this exploratory study, we employ a simple infrared dominant model for the gluon exchange 






that implements the vertex dressing. The resulting model, involving two distinct representations of 



f^ I the effective gluon exchange kernel, preserves both the axial-vector Ward-Takahashi identity and 

^^ ■ charge conjugation symmetry. Numerical results confirm that the pseudoscalar meson retains its 



Goldstone boson character. The vector meson mass, already at a very acceptable value at ladder 



Q^' level, receives only 30 MeV of attraction from this vertex dressing. For the axial-vector states, 

i-Q ' which are about 300 MeV too low in ladder approximation, the results are mixed: the 1"' state 

•'^ . receives 290 MeV of repulsion, but the 1^^ state is lowered further by 30 MeV. The exotic channels 

X, 

H , and 1 ^ are found to have no states below 1.5 GeV in this model. 
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I. INTRODUCTION 

In the modeling of QCD for hadron physics, the rainbow truncation of the quark prop- 
agator Dyson-Schwinger equation (DSE) coupled with the ladder truncation of the Bethe- 
Salpeter bound state equation (BSE) has been found to provide a very efficient descrip- 

riQ n 

tion for ground states ll|, |2j and for finite temperature and density p|. In particular, a 
ladder-rainbow model |^, with one infrared parameter to generate the empirically accept- 
able amount of dynamical chiral symmetry breaking, provides an excellent description of 
the ground state pseudoscalars and vectors including the charge form factors p, |6|, elec- 
troweak and strong decays J7|, |8| , and electroweak transitions J9|, |lO| • Recent comprehensive 
reviews of QCD analysis and modeling for nonperturbative physics emphasize the gauge 
sector DSEs |ll| and hadron physics ^ . The ladder-rainbow truncation is known to satisfy 
both the vector and the fiavor non-singlet axial-vector Ward-Takahashi identities; the lat- 
ter implementation of chiral symmetry guarantees the Goldstone boson nature of the fiavor 
non-singlet pseudoscalars independently of model details [l^. The small explicit symme- 
try breaking through current masses provides the detailed description of the pseudoscalar 
masses. 

Although the vector masses are not explicitly protected by a symmetry, the excellent 
description in ladder-rainbow truncation (typically within 5% of experiment i4|) illustrates 
the strong correlation between hyperfine splitting amongst "S-wave" states and dynami- 
cal chiral symmetry breaking. However, the ladder-rainbow truncation has inadequacies 
that are beginning to be understood and addressed. Corrections to the bare quark-gluon 
vertex inherent in the ladder-rainbow kernel have been examined within a schematic in- 

n 

frared dominant model 13] that admits algebraic analysis. In a dressed loop expansion, the 
corrections to the ladder-rainbow Bethe-Salpeter kernel were found to have repulsive and 
attractive terms that almost completely cancel for pseudoscalars and vectors but not so for 
scalars ll4\. Our understanding is far from complete; the schematic model used for this 
analysis does not bind scalar or axial vector meson states. 

In the axial vector channels 1^+ (01//1) and 1+^ (bi/hi), the ladder-rainbow models, 
exemplified by Refs. y,ll5|, are too attractive: the masses produced are 0.8-0.9 GeV [l6[. 
Evidently the orbital excitation energy rriai — rrip ^ is about a factor of 3 too small. On the 
other hand, very reasonable rua-^ and rrih-^ values (~ 1.3 GeV) are obtained with covariant 



separable models |l7l . llSl |l9| that incorporate some of the key features of the ladder-rainbow 
truncation. This encourages a systematic examination. 

The ladder Bethe-Salpeter kernel is vector-vector coupling (7^ ® "^u) and this generates 
a particular coupling of quark spin and orbital angular momentum. Some of the processes 
beyond this level constitute dressing of the quark-gluon vertex which generates a more 
general Dirac matrix structure for the vertex. Although 12 independent covariants are 
needed to describe the most general dressed vertex, those that involve the scalar Dirac 
matrix are of particular interest. Meson bound states are dominated by the infrared and one 
expects that it is most important to model the dressed vertex at very low gluon momentum. 
In this case only the quark momentum and the Dirac matrices are available for construction 
of the vector vertex covariants; of the three possible covariants, two involve 7^ and the 
third involves the Dirac scalar matrix. The latter generates a different coupling of quark 
spin and orbital angular momentum; such a correction to ladder-rainbow could distinguish 
between "P-wave" states such as the axial-vectors and the "S-wave" pseudoscalars and 
vectors. Another characteristic of a Dirac scalar matrix term in the vertex is that, in 
the chiral limit, it cannot be generated by any finite order of perturbation theory; it is 
generated by dynamical chiral symmetry breaking in quark propagators internal to the 
dressed vertex. Since dynamical chiral symmetry breaking sets the infrared scale of the 
ladder-rainbow truncation, it is natural to include the Dirac scalar part of the dressed 
vertex in considerations beyond this level. 

It is only recently that lattice-QCD has begun to provide information on the infrared 

structure of the dressed quark-gluon vertex [20;]. In the absence of well-motivated non- 

perturbative models for the vertex, many authors have employed the (Abelian) Ball-Chiu 

Ansatz 12 ll times the appropriate color matrix and comprehensive results from a trunca- 
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tion of the coupled gluon-ghost- quark DSEs have been obtained this way [22]. However, 
there is no known way to develop a BSE kernel that is dynamically matched to such quark 
propagator solutions in the sense that chiral symmetry is preserved through the axial- vector 
Ward-Takahashi identity. It is known that such a symmetry imposes a specific dynamical 



relation between the quark self-energy and the BSE kernel 23|. Since the ladder-rainbow 
kernel contains infrared phenomenology matched to the chiral condensate, a specification of 
correction terms that ignores chiral symmetry will needlessly alter the pseudoscalar sector. 
An artificial fine tuning of parameters can misrepresent the relationship to other sectors. 



There is a known constructive scheme [lj| that defines a diagrammatic expansion of the 
BSE kernel corresponding to any diagrammatic expansion of the quark self-energy such that 
the axial-vector Ward-Takahashi identity is preserved. We will apply this to a 1-loop model 
of the dressed vertex. Since use of a finite range effective gluon exchange kernel to construct 
the 2-loop BSE kernel leads to a very large computational task, especially with retention of 
all possible covariants, we begin here with a simplification. To an established finite range 
rainbow self-energy, we add a 1-loop vertex dressing model using the Munczek-Nemirovsky 
(MN) [l3| delta function model. The corresponding BSE kernel is then formed by the previ- 
ously mentioned constructive scheme, with one modification: to preserve charge conjugation 
symmetry for appropriate meson solutions, this approach requires that the BSE kernel be 
symmetrized with respect to interchange of the two distinct effective gluon exchange models 
that appear therein. The resulting kernel still preserves the flavor non-singlet axial-vector 
Ward-Takahashi identity. This hybrid model retains the advantages of a flnite range ladder- 
rainbow term, providing about 0.9 GeV towards axial-vector masses, while enabling a fea- 
sible exploration of the effects of vertex dressing. Recent investigations of the BSE beyond 



rainbow-ladder truncation |2J, |2^ have exploited the algebraic structure that follows from 
use of the MN model throughout and have not been able to address axial-vector states. 

In Section |n] we consider vertex dressing within the quark DSE and specify the employed 
rainbow self-energy and the model dressed vertex. The resulting 2-loop self-energy of the 
present approach is described. The chiral-symmetry-preserving BSE kernel is obtained from 
this self-energy in Section IIIII where the preservation of the axial-vector Ward-Takahashi 
identity in the present context is outlined. We discuss the solutions for the dressed quark 
propagator in Section HVl In Section|3we describe the BSE of the present work. Numerical 
results for meson masses are presented in Section IVTl Section FVIII contains a summary and 



the Appendix provides details of the 2-loop quark DSE that arises here. 

II. VERTEX DRESSING AND THE QUARK DYSON SCHWINGER EQUATION 

We work with the Euclidean metric wherein Hermitian Dirac 7-matrices obey {7^, 7,^} = 
2(5^j,, and scalar products of 4-vectors denote a ■ b = I]f=i ciibi. The color group is SU{Nc) 
with Nr = 3. 



In QCD, the Dyson-Schwinger equation for the renormahzed quark propagator is 

S-\p) = Z2{i^ + m) + ZiF J"^ dk g^D^^ip - k) ^If^Sik) yr,(A;,p) , (1) 

where dk = d'^k/{2'n')^, g is the renormahzed couphng constant, the A°/2 are the SU{3) 
color matrices, m is the quark bare mass, D^jj{q) is the renormahzed dressed gluon prop- 
agator in Landau gauge {q = p — k), and r,,{k,p) is the renormahzed dressed quark-gluon 
vertex. Here Zip and Z2 are the vertex and quark field renormalization constants. With 
a translationally invariant ultra-violet regularization of the integrals characterized by mass 
scale A, the renormalization conditions are S{p)'^ = i;^ + m^{fi) and T^{p,p) = •y^, at a. suf- 
ficiently large spacelike renormalization point p"^ = fi^. Here mR is the renormahzed current 
mass related to m through Z2m = Z/^rriii, with Z4 being the renormalization constant for 
the scalar component of the self-energy. 

Most studies have used the rainbow truncation in which the kernel of Eq. (^ becomes 

g^Df,^{q) ZiF yrj,(fc,p) -^ A^,^(g) ^l^y , (2) 

where A^,^(g) = t^^{q) A(g^), t^y{q) is the transverse projector, and A(g^) is an effective in- 
teraction. Due to chiral symmetry, there is a close dynamical connection between the kernel 



of Eq. (0) and the Bethe-Salpeter kernel for pseudoscalar mesons |l2l.ll4l|. This connection. 



and the observation that A(g^) should implement the leading renormalization group scaling 



of the gluon propagator, the quark propagator and the vertex, has been exploited [2y| to 
specify the ultraviolet behavior of A(g^) by that of the renormahzed quark-antiquark in- 
teraction or ladder Bethe-Salpeter kernel. In such a QCD renormalization group improved 
ladder-rainbow model J2^, A(g^) behaves as 47r Q;^^^°°P(g^)/g^ in the ultraviolet and is pa- 
rameterized in the infrared. The resulting rainbow Dyson-Schwinger equation reproduces 
the leading logarithmic behavior of the quark mass function in the perturbative spacelike 
region. The corresponding ladder Bethe-Salpeter kernel is 

\ a \a 

ir^-A^,(g)y7^ ®y7, . (3) 

Independently of the details of the model, this ladder-rainbow truncation preserves chiral 
symmetry as expressed in the axial-vector Ward-Takahashi identity, and thus guarantees 
massless pseudoscalars in the chiral limit [l^]- It is known that the chiral symmetry relation 



between the kernels of the Dyson-Schwinger and Bethe-Salpeter equations may be main- 
tained order by order beyond the ladder-rainbow level in a constructive scheme in which the 
first two terms are |l4| 

Zip ^T,{k, p) = ^^^ - J dl A,,(/) ^7p S{k + /) ^ 7. S{p + I) y7A • (4) 

This replaces the corresponding portion of the Dyson-Schwinger kernel of Eq. 0. Normally 
one has A = A. 

The corresponding Bethe-Salpeter kernel in this scheme is generated as the sum of terms 
produced by cutting a quark line in the diagrammatic representation of the quark self-energy. 
Due to the complexity of a Bethe-Salpeter calculation with a two-loop kernel, existing nu- 



merical implementations 1J,|2J,|2^ for bound states have used the infrared dominant model 
in which A(g^) = A(g^) oc S^{q). Through its results for dressed quark propagators and 
pseudoscalar and vector mesons, it is known that this simple model effectively summarizes 
the qualitative behavior of more realistic models. The strong infrared enhancement is the 
dominant common feature. The algebraic structure that this model gives to the DSE-BSE 
system was exploited to implement the BSE kernel consistent with an Abelian-like summed 
ladder model for the quark-gluon vertex |2^. It was found that the quark propagator and 
the BSE meson masses are well-represented by the 1-loop version of that vertex model. A 
disadvantage of this algebraic model for the entire BSE kernel is that it does not support 
non-zero relative momentum for quark and anti-quark. This is likely the reason why it does 
not bind the (P-wave) axial-vector states. 

Since existing finite range ladder-rainbow models typically produce axial vector masses of 
about 0.9 GeV, we use this to build a model that generalizes the approach of Ref. [1^; the 
delta function interaction is used only as an effective representation of the gluon exchange 
that implements vertex dressing via Eq. (j3)). Since this hybrid model has A 7^ A, an adjust- 
ment must be made in the construction of the symmetry-preserving Bethe-Salpeter kernel 
and we address this in Section UTTl 

Hence for A^^{q) appropriate to the ladder-rainbow level of the present approach, we 



take the following convenient representation 



m 



A^,(g)=47r2DV(g)^exp -^ . (5) 



The parameters D and u are required to fit the light pseudoscalar meson data and the 
chiral condensate; throughout we will use the values |l5| tu = 0.5 GeV, D = 16 GeV^^. We 








FIG. 1: The approximation employed here for the dressed quark-gluon vertex. To provide for 
manageable calculations, the internal gluon line that implements the dressing is simplified to the 
5- function form A given in Eq. ©, and identified in diagrams as a wavy line. 

shall consider only w/fi-quarks with m = 5 MeV representing the current mass Il5|- With 

nnn 

this form we retain the phenomenological successes of recent studies j^, |15|, |26|] of the light 
quark flavor non-singlet mesons at the rainbow-ladder level. We note that the above model 
does not implement the ultraviolet behavior of the QCD running coupling; this contributes 
typically 10% or less to meson masses |4] and this level of precision is not our concern here. 
The form of this model interaction also produces ultraviolet convergent integrals and thus 
renormalization will be unnecessary here; one has Zip = Z2 = 1 and m^ = m. 

For the effective gluon exchange that implements vertex dressing, we employ the Munczek- 
Nemirovsky (MN) model [l3| in the form 

Apx{l) = (271)^ G6\l)t,x{l) . (6) 

The single parameter G represents the integrated infrared strength which is the key feature 
for empirically successful ladder-rainbow models. To set G, one could match the quark mass 
function at p^ = or the vector meson ground state mass in rainbow-ladder truncation. Here 
we have chosen to simply demand that both models have the same integrated strength. From 
Eqs. © and (jHI), this gives G = Duj^/2, and yields G = 0.5 GeV^. We note that the smaller 
estimate G ~ 0.25 GeV^ is suggested by reproduction of mp = 0.770 GeV in rainbow-ladder 
truncation. We present our results for the parameter range G = — 0.5 GeV^. 

The simple form, Eq. (jH)), for A reduces Eq. (j3]) for the dressed vertex to an algebraic 
form. With X^ /2X'' /2\^ /2 = (-1/6) AV2, one obtains 

G 

^u{k,p)=iy + —ipS{k)-i^S{p)-fp , (7) 

where the factor 1/8 comes from the above 1/6 times an extra factor of 3/4 generated by 
the combination of the transverse projector and the 5 function. With this, Eq. (Q) for the 



quark propagator becomes 

S-\p) = ^^ + m + J dkA,,{p-k)i^^^,S{k)^. + ^^,S{k)^,S{k)^,S{p)^pj . (8) 

Our approximation for the dressed quark-gluon vertex is summarized by Fig.^where the 
wavy hue identifies the gluon exchange that has been simphfied to the MN model. We note 
that, although in an Abelian theory like QED this is the only 1-loop diagram that provides 
vertex dressing, this is not the case in QCD where there is a second 1-loop quark-gluon 
diagram allowed due the existence of a 3-gluon vertex. There is no definitive information 
available for nonperturbative modeling of the 3-gluon vertex and such considerations are 
beyond the scope of this work. We note that in Ref. 25| the effect of such a contribution 
was estimated by arguing that a rescaling of the Abelian-like 1-loop diagram for the quark- 
gluon vertex would be the result. The examination of Bethe-Salpeter bound states in this 
context again utilized the algebraic structure afforded by the MN model. We explore the 
present model for its capability to addresss a larger class of meson states (e.g., axial- vectors) 
through an extension of Ref. jlj] where the vertex phenomenology is based upon Eq. Q. 
The establishment of this capability may allow a wider examination of vertex phenomenology 
in future investigations. 

III. SYMMETRY-PRESERVING BETHE SALPETER KERNEL 

The close dynamical connection between the Bethe-Salpeter kernel and the quark propa- 
gator is manifest in chiral symmetry as expressed through the axial- vector Ward-Takahashi 
identity. We shall be interested only in the color singlet, flavor non-singlet channels where 
such an identity leads to the Goldstone phenomenon in the chiral limit. The flavor singlet 
channels have an axial anomaly term in the axial-vector Ward-Takahashi identity, which 
blocks the Goldstone phenomenon; the t] — rj' system and scalars with flavor singlet compo- 
nents lie outside of our considerations. From the quark self-energy of our model as given in 
Eq. (jHl), the symmetry-preserving BSE kernel can be obtained by the constructive scheme 
of Ref. Il4| with a generalization to account for the concurrent use of two distinct effective 
interactions. 

In a flavor non-singlet channel, and with equal mass quarks, the axial-vector Ward- 



Takahashi identity is 

-iP^Tl{p;P) = S-\p+h, + ^,S-\p^)-2m^T'{p;P) , (9) 

where we have factored out the exphcit flavor matrix. The color-singlet quantities F^ and F^ 
are the axial- vector vertex and the pseudoscalar vertex respectively, p is the relative quark- 
antiquark momentum, and P is the total momentum. We use the notation p+ = p + C,P and 
p^ = p— (1 — ^)P, where the momentum sharing parameter ^ = [0, 1] represents the freedom 
of choice of relative momentum in terms of individual momenta. Due to Poincare invariance, 
results will be independent of ^ if a complete momentum and Dirac matrix representation 
is used J5|. In the present work with equal mass quarks, we use ^ = 1/2 for convenience. 

The BSE kernel determines the dressed parts of the vertices F|^ and F^, while the self- 
energy determines the dressed part of S^^. Corresponding to a given approximation for S~^, 
one seeks the matching approximation for the BSE kernel so that the above Ward- Takahashi 
identity holds for the truncated theory. In this way the pseudoscalar bound state results 
will be dictated by the pattern of chiral symmetry breaking and largely invariant to model 
details. 

If we deflne 

A{p;P) = -zP^Tl{p-P) + 2m^r'{p-P) , (10) 

then the BSE integral equations for F^ and F^, that have inhomogeneous terms 2^2 757/^ and 
.^4 75 respectively, can be combined to yield 

A(p; P)„^ = [~tZ2 75f + ^4 2mR75]«^ + j'' dkK{p, k; P),^,,^[5(A:+)A(A;; P)S{k_)]^s , 

(11) 
where K{p, k; P)af3,Sj is the BSE kernel. For generality and clarity we have retained the 
renormalization constants although, in the eventual application to the present model, they 
will be unity. If Eq. Q is to hold, then we may eliminate A{p; P) in favor of S^^ so that 
Eq. ()11|1 explicitly relates the BSE kernel and the dressed quark propagator via 



[S ^(p+)75 + 75S' ^{p-)]a/3 = ^2[-n5f + 2m75]„/3 

+ 1"^ dkK{p,k;P)^p,s-yb,S{k_)+S{k+h,]-,s . (12) 

It is helpful to move the inhomogeneous term to the left hand side so that cancellations leave 
only the self-energy integrals on the left. With use of the explicit form of the self-energy 



integrals from the DSE, as given by Eq. (^, we thus have the axial- vector Ward-Takahashi 
identity expressed in the equivalent form 

Zw- J dk\^g^D^^{p+ - k)['yf,S{k)T^{k,p+)-f5]^i3 + g^D^^{p_ - k)['y5jf,S{k)T^{k,p_)]af3\ 

= p dkK{p, k; PU,s-y[l5S{k^) + S{k+h,]^s . (13) 

For a given approximation or truncation to the vertex r,^(fc,p), the corresponding truncated 
BSE kernel has to satisfy this integral relation to preserve chiral symmetry. Any proposed 
Ansatz can be checked through substitution. The rainbow truncation on the left (substitute 
Eq. 0) and the ladder truncation on the right (substitute Eq. Q) obviously satisfy Eq. (fT!^ . 
The general relation between the BSE kernel K and the quark self-energy E can be 



expressed through the functional derivative 



23| 



K{x\y';x,y) = -—^j:ix',y') . (14) 

dS{x,y) 

It is to be understood that this procedure is defined in the presence of a bilocal external 
source for qq and thus S and S are not translationally invariant until the source is set to 
zero after the differentiation. An appropriate formulation is the Cornwall- Jackiw-Tomboulis 

n 

effective action [22l|. In this context, the above coordinate space formulation ensures the 
correct number of independent space-time variables will be manifest. Fourier transformation 
of that 4-point function to momentum representation produces K{p, q; P) having the correct 
momentum flow appropriate to theBSE kernel for total momentum P. 

The constructive scheme of Ref . [ij] is an example of this relation as applied order by order 
to a Feynman diagram expansion for S(p). An internal quark propagator S{q) is removed 
and the momentum flow is adjusted to account for injection of momentum P at that point. 
With a change in sign (related to use of {75,7^} = in Eq. (fT!?jl ). this provides one term 
of the BSE kernel. The number of such contributions coming from one self-energy diagram 
is the number of internal quark propagators. Hence the rainbow self-energy generates the 
ladder BSE kernel; this is the first term in Fig. |21 A 2- loop self-energy diagram (i.e., from 
1-loop vertex dressing as in the present case) generates 3 terms for the BSE kernel. By 
substitution into Eq. ()1H|1 one can confirm that the axial-vector Ward-Takahashi identity is 
preserved. Similarly, the vector Ward-Takahashi identity is also preserved. 

In the present study, our combined use of two different gluon propagator models requires 
that the BSE kernel obtained so far be subjected to an additional procedure to preserve 

10 



charge conjugation symmetry. The behavior of the Dyson-Schwinger equation under charge 
conjugation leads to invariance of the self-energy amplitudes under interchange of the two 
different gluon propagators. In particular, Eq. (jSj) may be written in the equivalent form 

S-\p)=^^ + m + J dkA^,{p-k)i^^^,S{k)^, + j^,S{p)^,S{k)^^S{k)^,'j , (15) 

and it makes no difference which of the two gluon propagator models is considered to be 

implementing the 1-loop vertex dressing. However, once a quark propagator is removed 

from a self-energy diagram to produce a contribution to the BSE kernel, interchange of the 

two gluon lines, or reversal of the quark line directions, produces a distinct contribution 

to the BSE kernel. The required invariance of the BSE kernel to charge conjugation can 

be restored if the 3 diagrams are symmetrized with respect to interchange of gluon lines, a 

procedure that does not alter the quark propagator. Thus with the present hybrid model 

there are three pairs of 2-loop diagrams to be used for the BSE kernel; these are displayed 

in Fig. H 

With these considerations, the explicit form for the two-loop BSE kernel of the present 

model is 

4 
K{p, k; P)ap;5j = -g A^^(p - k) [7^]^^ [7^]^^ 

-g J dq [A^^(p - k)Apxiq) + A^^(p - k)Apxiq)] b^,]^^ [ipSik. + q)-f^S{p_ + qhxlg^ 

-g J dq [A^^(p - k)Apx{q) + A^j,(p - k)Apx{q)\ bpS{p+ + q)-fpS{k+ + q)-fx]^^ [-fu]sp 

-g y rfg [Apx{q)A^u{k ~p-q)+ Apx{q)Ap^{k - p - q) x 

bpS{p+ + q)lp]^^ b\S{k_ - q)lu]^p . (16) 

It is straightforward to verify by direct substitution of this kernel and the 1-loop dressed 
vertex that the axial-vector Ward-Takahashi as expressed by Eq. (fTT^ is satisfied. The charge 
conjugation invariance of the kernel can be verified by substitution into the bound state BSE 

T{p;P)c,p = J dkK{p,k]P)ap.STX{k]P)-yS , (17) 

followed by charge conjugation in the form 

CT^{-p;P)C-'=CMr{p;P) . (18) 

Here C = —7274 is the charge conjugation matrix. Cm is the charge parity of the meson, 
and x(p; P) = 'S'(p_|_) T{p; P) S{p_) is the Bethe-Salpeter wave function. 

11 
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FIG. 2: The 2-loop Bethe-Salpeter kernel of the hybrid model that utilizes two different represen- 
tations of effective gluon exchange: A denoted by springs, and A denoted by wavy lines. 

IV. VERTEX CORRECTION TO THE QUARK PROPAGATOR 

The dressed quark propagator can be represented by a pair of amplitudes and the two 
convenient forms used herein are defined by 

S{p) = [zMip') + B{p')]-' = -ii>ay{v^) + as{v^) . (19) 

Projection of the Dyson-Schwinger equation (jH)) of the present model on to the amplitudes 
A and B yields the coupled equations 

where x = j? . In Appendix 1X1 explicit expressions are given for liix) and Jj(x) as integrals 
involving oyiy) and asiy) for all spacelike y = k"^. For z = 2,3, the integrals /j, Jj are 
proportional to the strength G for vertex dressing. In the above they multiply scalar and 
vector propagator amplitudes at x because, due to the delta function, one quark propagator 
involved in the 2-loop self-energy is always evaluated at the external momentum. The quark 
Dyson-Schwinger equation thus has partly an algebraic structure and partly an integral 
structure. 
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The rainbow truncation, G ^ 0, yields A{x) = Ii{x) and B{x) = Ji{x); the nonhnearity 
here is contained in the expressions for /i(x) and Ji{x) given in Appendix ^ In this hmit, 
our results reduce to those of Ref . [l5[ thus providing a useful check of our solutions. The 
limit uj^ ^ provides a complementary check: the distinct effective interactions of the 
present model, A^^{q) and A^j,(g), become identical. In this case the quark DSE, Eqs. (f^ 
and ()21|) . takes on the algebraic structure that follows from use of the MN delta-function 
model in all aspects; this has been applied in Ref. [1^ and in the 1-loop vertex considerations 
of Ref. y. 

The spacelike (x > 0) solutions of Eqs. (j^ and PT|) are calculated by iteration subject 
to the boundary conditions A{x) -^ 1 and B{x) —^ m for large enough spacelike x. We note 
that the solution method can be based on purely numerical iteration or upon a mixture of 
numerical iteration followed by determination of polynomial roots. The latter case empha- 
sizes the algebraic dependence of Eqs. (^Uj) and (PT|) upon the explicit A{x) and B{x) once 
the integrals /« and Jj have been evaluated. Both methods were used as a check. 

The amplitudes A, B in the spacelike region are shown in Fig. IHl for the chiral (m = 0) 
limit. The rainbow result is displayed as a solid line and is compared to results that include 
vertex dressing characterized by strength parameter G, with G = 0.5 GeV^ being our esti- 
mate of the internally consistent value for this model. The influence of vertex dressing is 
evident mainly in the infrared region x < 1 GeV^ where both A{x) and B are reduced with- 
out an alteration in the qualitative behavior. There is a modest increase in the dynamical 
mass B/A. These results parallel those obtained from the MN delta-function model J2J|. 
At X = 0, B decreases by ~ 20% and A — 1 decreases by ~ 30%. 

The effect of vertex dressing on the quark propagator can be cast in a more physical con- 
text by considering the chiral vacuum quark condensate (qq)^ which characterizes dynamical 
chiral symmetry breaking. In general the condensate at renormalization scale fi is 

A 



(gg)° = - hm Z,{fi',A^) N.Tm f dkSo{k) , (22) 

;he chiral limit propagator. For an extensive discussion of the condensate in 
2S|. The present simple model is ultraviolet convergent, regularization is not 



where 5*0 (A;) is 
QCD, see Ref. 
required, and Z4 = 1. Eq. ([2^ then yields 



„.»o_ 3 r,.... B„(y) 



»«) =-4^/„ ^'\A„(yY-,B,(yY ' <"'' 
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FIG. 3: Spacelike quark propagator functions A(x) (top curves) and B{x) (lower curves) in the 
chiral limit, m, = 0. The rainbow truncation (solid line) is compared to the 1-loop dressed vertex 
model with two representative strength parameters. 

which can be considered to be at the typical hadronic scale /i ~ 1 GeV. The results for 
(gg)° given in Table IJ show a decrease of 2% as the rainbow truncation (G = 0.0 GeV^) is 
extended by the range of vertex dressing considered here. 

Since the 2-loop BSE kernel used here preserves the axial- vector Ward-Takahashi identity, 

I I 

the resulting mass relation [12] for pseudoscalar mesons is also preserved. At low ruq this 
relation becomes the GMOR relation 



{f:ymi = -2m,{qq)' + 0im: 



(24) 



where /° = /°/v2, and /^ is the chiral limit leptonic decay constant in the convention 
where Z^^^* = 0.131 GeV at the physical u/d-qnaik mass. From previous explorations [IJ, 
|2J], we expect m^r to be quite stable to vertex dressing; the relative insensitivity of (gg)° 
should therefore be evident in /^ also. Here we estimate f2 by using the dominant term 
i'y^B(j{q'^)/ f^ of the chiral pion Bethe-Salpeter amplitude J29| 



u. 



0\2 



2aJ dqBoiq' 



2772 70 



dyyBo{y) 



((t',(Ts — cr^cr^) + a^CTs 



{<yWs - (^v<^'s) + (^v<ys 



(25) 



This expression is known to provide an underestimate by about 10%. 
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The results for /^ are also given in Table HI The variation between the 1-loop (rainbow) 
kernel {G = 0.0 GeV^) and the estimate for the physical 2-loop kernel (G = 0.5 GeV^) is 
1%. The expected stability of m^ is confirmed from the BSE solution in Sec. El 

TABLE I: The chiral quark condensate (qq)^ and the chiral leptonic decay constant /^ for a range 
of vertex dressing strength G (in GeV^). 



G 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


Expt 


(_(^g)0)l/3 


0.2511 


0.2505 


0.2498 


0.2492 


0.2485 


0.2478 


0.22-0.24 GeV 


/° 


0.1190 


0.1188 


0.1185 


0.1183 


0.1180 


0.1178 


0.131 GeV 



To facilitate our investigation of the BSE with the 2-loop kernel, the continuation to 
P^ = — M^ for total meson momentum will be implemented under the real axis approx- 
imation [7], defined as the substitution F{x) —>■ F{^x) where F{x) stands for the quark 
propagator amplitudes A{x) and B{x) as they appear in the BSE kernel. The otherwise 
necessary, model-exact, method requires knowledge of the complex momentum plane struc- 
ture of the BSE and DSE kernels, and has been developed and tested |y| only for ladder 
truncation. Since the present 2-loop kernel is exploratory and future developments are ex- 
pected, we use the approximate method here. With a ladder kernel of the present type, the 
real axis approximation produces values for rup, itik* and ms that are within 4%, 7% and 

n 

1% respectively f?] of the model-exact values; an accuracy of about 10% is anticipated for 
masses up to about 2 GeV J7|. 

Propagator amplitudes A and i?, in the chiral limit, and along the real timelike p^ axis 
are displayed in Fig. E] for representative values of the vertex dressing strength G. For equal 
mass quarks, the meson mass M and the quark x = p^ required by the BSE kernel along the 
real axis are related by x > — M^/4. Thus the timelike range shown in Fig. |3]is sufficient 
for a meson mass up to 2 GeV. The repulsive effect of our model vertex dressing is evident 
in the enhanced dynamical chiral symmetry breaking [B > 0) in the timelike region. In 
models of this type, the strong dynamical chiral symmetry breaking creates a significant 
infrared timelike domain where the dressed quark has no real mass shell; qq bound states 
with masses within this domain have no spurious qq decay mode. We illustrate this in Fig.E] 
which displays the propagator denominator xA{x)'^ + B{xy in the chiral limit. A zero of this 
quantity indicates a quark mass shell. In rainbow truncation (G = 0.0), the physical domain 

15 




3 

_ 2 

> 

O 1 

CN 

+ 

CN 

< 

X 

-1 

-2 



, 1 1 1 1 1 1 1 

\ 


r* 


.^-" — ~-^--i;:i£s-'^'*^ 


- - G=O.OGeV^ 
G=0.01 GeV^ 

G=0.1GeV^ 

G=0.3 GeV" 

G=0.5GeV" 


- 




' , 1 


1 




1 -0.5 


0.5 


1 



X (GeV ) 



X (GeV ) 



FIG. 4: The influence of vertex dressing {G > 0) on the chiral hmit quark propagator for timehke 
and spacehke real x = p^. Left Panel: Propagator amphtudes A{x) and B{x). The sohd hne is 
a good representation of the rainbow truncation. Right Panel: The denominator function. Finite 
vertex dressing strength prevents a zero (a physical mass-shell) in this timelike region. 

of applicability identified this way is x > —0.7 GeV^ [15], which allows meson states below 
about 1.7 GeV to be free of spurious qq widths. In Fig. HI our model vertex dressing can be 
seen to enlarge the physical domain of applicability; any meson states up to 2 GeV would 
be without such spurious widths. This is consistent with the effective increase in infrared 
strength that is generated by vertex dressing. 

V. THE BETHE SALPETER EQUATION AND MESONS 

With the constructed 2- loop kernel, Eq. (fTB|). the homogeneous Bethe-Salpeter equation 
for mesons becomes 



T{p;P) = / d k K{p, k; P)ap-5'yXik; P)'y5 
4 f 



G 

12 
G_ 

12 

Q 

— I dkA^^{p- k) {7px(p; P)li,S{k^)-ipS{k_)-i^ + 7^S(A;+)7p5(A;+)7j.x(p; P)lp} 



dkAf,^{p - k) {-ff,xik; P)-fpS{k^)-fuS{p-)-fp + -fpS{p+)-fpS{k+)jpxik; P)-fu} 
dkApy{p - k) {-ipS{k+)-ipx{k] P)-f^S{p^)-fp + -fpS{p+)-fpx{k; PhpS{k^)-f^} 



(26) 
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where the integrations over the (5-functions have been carried out and xiP] P) = 
S'(p+)r(p; P)S{p-). Discrete solutions for meson mass M„ exist for P^ = — M^. 

To solve for a particular meson, one must specify the appropriate quantum numbers J^*-^ 
by expressing the bound state vertex F in the general covariant form that has the corre- 
sponding transformation properties. The quark flavor content can be specified by the current 
quark masses. Here we consider pseudoscalar, vector and axial- vector charge eigenstates and 
adopt isospin symmetry through degenerate u and d quarks. The general covariant forms 
used here are 

r^^(p; P) = 75 [ro(p; P) - ifTiip; P) - ifV^^v. P) - [f,f\ TsIp; P)] , (27) 

r;:(p; P) = il bro(p; p) + f ri(p; p) - i>V2{p; p) + 1 W.i>] TsIp; p)] 

+pl[T2{p-P) + 2irT^{p-P)] 



+pI [T,{p- P) + ifT,{p- P) - ii)T,{p- P) + [f,i}] Ty{p- P)] , (28) 

T^\p;P) = ^,jj:[tToip;P) + r^iip;P)-^T,ip;P) + t[rj]Tsip;P)] 
+l5pl[r2{p;P) + 2zrr3{p;P)] 
+15pI [Ta{p; P) + ^fT^{p- P) - ti>T,{p- P) + [f,;^] Tj{p- P)] , (29) 

where the notation a^ denotes a 4-vector transverse to P^, i.e., a^ = t^y{P)ay. The vector 
and axial-vector amplitudes Tj^/"^^ are transverse to the total momentum P. The Fj are 
scalar functions of p^, p-P. For charge eigenstates, the Fj are either odd or even under the 
interchange p-P ^ —p-P. 

For numerical solution of Eq. ^I^ we project onto the basis of Dirac covariants to obtain 
coupled equations for the amplitudes Fj. The dimensionality of integration is reduced by 
expansion of Fj in the complete orthonormal set of Chebyshev polynomials Tk{z) in the 
variable z = p-P/y/jP^, 

F,(p^ z- P') = E ^' T,{z) Tlip'; P') . (30) 

fc=0 

Projection of the equation onto the T^ basis then produces a set of coupled 1-dimensional 
integral equations for the F^(p^; P^). An explicit factor of p-P is extracted from the odd 
Fj and an even order Chebyshev expansion is used for the remaining factor. Convergence 
and stability is normally achieved with Nch = 4 or 6. The meson mass M„ is determined by 
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introducing a linear eigenvalue A(M) so that the Bethe-Salpeter equation reads 



A(M)r(p;P) 



a/3 



dk K^p,s-y{p,k-P) [S{k.)T{k;P)S{k^)] 



■yS 



(31) 



with P^ = — M^; the mass M is varied until X{M = Mn) = 1. The physical ground state in 
a given channel is determined by the largest real eigenvalue. 

VI. NUMERICAL RESULTS AND DISCUSSION 

We explore ground state u/d-qnaik mesons and adopt the ladder kernel model of Ref. [15| 
which provides the parameters u = 0.5 GeV, D = 16 GeV~^ (see Eq. (0)) and rriu/d = 5 MeV. 
The strength parameter G for vertex dressing is varied in the range 0.0-0.5 GeV^. Our 
model Bethe-Salpeter kernel produces real ground state masses since it does not contain a 
mechanism for a hadronic decay width; it also does not distinguish between isoscalar and 
isovector states. We indicate the physical state identifications with these qualifications in 
mind. Real mass solutions of the type obtained here have in the past provided a basis 
for successful perturbative descriptions of strong decay widths in non-scalar channels, an 
example being the p — > tttt width j^, l3lLl82l|. 

TABLE II: The 0~+ (vr) and 1"" (p) masses in GeV. 



G (GeV2) 





0.1 


0.2 


0.3 


0.4 


0.5 


Expt 


iVch = 2 m„ 


0.143 


0.142 


0.142 


0.141 


0.141 


0.140 




nip 


0.806 


0.808 


0.810 


0.811 


0.812 


0.813 




iVch = 4 m^ 


0.140 


0.140 


0.140 


0.140 


0.140 


0.139 




rup 


0.794 


0.791 


0.788 


0.783 


0.777 


0.771 




iVch = 6 m^ 


0.140 


0.140 


0.140 


0.140 


0.141 


0.139 


0.139 


nip 


0.794 


0.790 


0.785 


0.779 


0.772 


0.763 


0.770 



Table El records the pion (0~^) mass values for vertex dressing strengths between 
G = 0.0 GeV^ (rainbow-ladder truncation) and G = 0.5 GeV^ (our estimate of the physical 
value). The mass is practically insensitive to model details; this illustration of chiral sym- 
metry preservation through the axial-vector Ward-Takahashi identity has been observed in 
related earlier work with a simplified ladder kernel IJ, |2J, |25| . Table |n] also illustrates the 
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stability and convergence with respect to A'^ch (the number of Chebyschev polynomials used 
to expand the angle dependence of the Bethe-Salpeter amplitudes according to Eq. (jHUj) ). 
This and other numerical techniques were applied until an accuracy of better than 5% was 
achieved for the masses. 
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FIG. 5: The rUg dependence of the pion mass near the chiral limit. The fitted dashed line confirms 
that the dressed vertex model (characterized by G = 0.5 GeV^) satisfies the GMOR relation, 
Eq. dHJ. 

With strength (7 = 0.5 GeV^ for vertex dressing, the mg-dependence of the pseudoscalar 
bound state mass near the chiral limit is displayed in Fig. |S| The form 



m. 



ai^TUg + 0-2 TT^q 



(32) 



provides a good fit with ai = 1.94 GeV^, and a2 = 0.231. The dominance of the first 
term (by a factor of about 10^) suggests that the GMOR mass relation, Eq. ()24|). has been 
reproduced. This is reinforced by noting that, according to Eq. ()24|1 . we should have the 
correspondance ai ~ 2w|(gg)°|//° and with the model value of (gg)° presented in Table HI 
we deduce from the fitted ai that /° = 0.1272 GeV. This chiral value is 3% lower than the 
experimental value /^ = 0.1307 GeV (at the physical mass) and is consistent with previous 
findings |26|. It is also consistent with the independent calculation of /° in Section HVl 

Fig. ini displays the eigenvalue for the exotic channel where the ladder kernel produces 
a mass of 1.38 GeV. Our model vertex dressing has a large repulsive effect; the physical value 



19 



c< 



a 



1 . 1 1 1 1 1 . 


1 ' 1 ' 




0" 

N =6 

ch 


— G=0 
G^O.lGeV^ 

— G=0.3 GeV^ 

— -- G=0.5 GeV^ 


- 




^ 


/ 


1 1 1 1 1 1 1 1 1 1 1 1 1 



0.2 0.4 0.6 0.8 

M 



1.2 1.4 



FIG. 6: Exotic pseudoscalar (0 ) meson Bethe-Salpeter eigenvalue from the ladder kernel {G = 0) 
and from the dressed vertex model (G > 0). 

for (7 = 0.5 GeV suggests that any mass solution would be above 2 GeV. There is a scalar 0"'""'" 
solution to the present model at a mass of 0.7 GeV in ladder approximation; this decreases 
to 0.6 GeV for the physical 2-loop kernel. The ladder value is consistent with previous 
models of this type |l^ llSl |22| • A physical identification for the 0"*"^ solution here is not 
appropriate because hadronic configurations such as KK are expected to play an important 
role, and the ladder-rainbow truncation is known to be deficient in this channel |l4| . 

The mass values for the 1 (p/^) channel are presented in Table ITTl for the relevant range 
of parameters G and Net- The effect of the present model of 1-loop vertex dressing is seen 
to be an attraction of only 30 MeV. In contrast, a repulsive change of about 70 MeV applies 
when the delta function representation is employed for all aspects of the 2-loop kernel J2J| . 
That simplification also provided an extension to the complete ladder summation of the 
one-gluon exchange mechanism for the vertex; the result being just a 75 MeV repulsion 24 1. 



Evidently, finite momentum range effects in the kernel can influence the net result from 
vertex dressing and whether there is net attraction or repulsion. 

However, the larger issue in the present context is why the typical effect beyond rainbow- 
ladder seems to be no more than ±10% for the ground state vector mass when there is 
no explicit protection from an underlying symmetry. With the pseudoscalar state fixed by 
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FIG. 7: Axialvector meson Bethe-Salpeter eigenvalue from the ladder kernel {G = 0) and from the 
dressed vertex model (G > 0). Left Panel: The l"*""*" (ai//i) state. Right Panel: The 1"' (bi/hi) 
state. 

chiral symmetry, the hyperfine splitting that can be generated by modeUng the quark-gluon 
vertex by effective gluon exchange is evidently a weak perturbation on the ladder kernel. 
Recent work indicates that significant attraction is provided to the dressed gluon-quark 
vertex by the triple-gluon vertex, and a schematic model implementation suggests tha t up 
to a 30% attractive effect on rup can be provided in this way beyond rainbow-ladder 25|. 
Further work is required for an understanding of how best to model the quark-gluon vertex 
for hadron physics. 

No mass solution below 2 GeV was found in the exotic channel l""*". In this mass range, 
the largest eigenvalue had reached 0.5-0.8. Masses above 2 GeV might be possible in a 
model of this type. 

The eigenvalue behavior for the axial vector solutions in the 1"'""'" (ai//i) and 1^^ {bi/hi) 
channels are displayed in Fig.[7|and also in Table ITTTl In previous work, the ladder truncation, 
constrained by chiral data, is generally found to be 200-400 MeV too attractive for these P- 



wave states [7, 



15|,|l6|. Our present results agree with this. The l"*""^ channel shows a 30 MeV 



of attraction due to the effect of 1-loop dressing added to the ladder kernel. However, in the 
1+^ {hi/hi) channel, we find a repulsive effect of 290 MeV above the ladder kernel result, 
yielding a value close to m^^ . We are not able to compare these findings with previous work 
on vertex dressing since such P-wave states do not have solutions in the models considered 
p..eviou.>v fo. tKat purpose Q Q Q. Oth. .tud.es of . a.d , based on the >adde. 
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TABLE III: The 1++ (ai//i) and 1+- (61 //ii) masses in GeV. 



G (GeV2) 





0.1 


0.2 


0.3 


0.4 


0.5 


Expt 


iVch = 2 nia^ 


1.061 


1.062 


1.060 


1.056 


1.050 


1.043 




rrib^ 


0.968 


0.983 


0.999 


1.017 


1.038 


1.064 




Nch = 4 rua^ 


0.989 


0.985 


0.980 


0.975 


0.969 


0.964 




rub^ 


0.952 


0.978 


1.008 


1.046 


1.094 


1.155 




iVch = 6 m„, 


0.989 


0.982 


0.975 


0.969 


0.963 


0.959 


1.230 


nib^ 


0.954 


0.983 


1.020 


1.070 


1.138 


1.227 




Nch = 8 TUb^ 


0.954 


0.984 


1.024 


1.080 


1.161 


1.243 


1.230 



rainbow truncation have used a separable 



the phenomenological instruments 



13 



la 



approximation where the quark propagators are 



l^ . these studies find more acceptable masses 



for both states in the vicinity of 1.3 GeV. 

The typical widths of the ground state axial vector mesons are about 20% of the mass 
and Svr and Air decay channels are prominent. Although widths of this magnitude have 
been successfully generated perturbatively for vectors 7] and axial-vectors il9] from BSE 
solutions that do not have the decay channels within the kernel, it is possible that such 
dynamics are responsible for important contributions to the masses beyond ladder-rainbow 
truncation. 

VII. SUMMARY 

The effect of a 1-loop model of quark-gluon vertex dressing on the masses of the light 
quark pseudoscalar, vector and axial-vector mesons has been studied. To facilitate cal- 
culations the vertex model consists of an effective single gluon exchange represented by a 
momentum 5-function with strength G. The ladder-rainbow term is represented by a fi- 
nite range effective gluon exchange. The axial-vector Ward-Takahashi identity and charge 
conjugation properties are used to construct a consistent Bethe-Salpeter kernel from the 
quark self-energy. This model extends recent explorations beyond the rainbow-ladder level 
by allowing axial-vector meson solutions. 

At spacelike momenta the vertex dressing results in only slight corrections to the quark 
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propagator. The Bethe-Salpeter equation was solved numerically for the ground state pseu- 
doscalar, vector and axial-vector mesons. The results show that corrections to the rainbow- 
ladder truncation generate very little change to pseudoscalar masses (and effectively vector 
masses) when they are defined consistently with chiral symmetry. 

The axial- vector mesons respond to the present vertex dressing model in a way that calls 
for further study. The mass of the 1+^ {hi/ hi) state, too small by about 300 MeV in rainbow- 
ladder truncation, was raised by 290 MeV by the vertex dressing thus giving a satisfactory 
value. However the 1++ (ai//i) state decreased in mass by about 30 MeV, leaving it still 
about 300 MeV too light compared to experiment. In the exotic meson channels, and 
1^^, the vertex dressing effects are repulsive; possible solutions are well above 1.5 GeV and 
beyond the range of the methods of the present study. 

There is very little information on the non-perturbative structure of the dressed quark- 
gluon vertex that can be used to guide a practical phenomenology. The Ball-Chiu (Abelian) 
Ansatz [2l|, often-used to generate the quark self-energy, is not useful here because the 
only known way to define a chiral-symmetry-preserving BSE kernel requires an explicit 
Feynman diagram representation of the self-energy. The present model adopts a one gluon 
exchange vertex structure motivated by a previous study. Dynamical information on the 
non-perturbative structure of the triple gluon vertex, and the contribution it can make to the 
quark-gluon vertex, is needed to further clarify how best to model the BSE kernel beyond 
ladder-rainbow. The attractive influence J2a] of the triple gluon vertex can have impor- 
tant consequences for modeling of the hadron spectrum. The high mass meson states are 
furthermore exptected to move significantly when including their hadronic decay channels. 

APPENDIX A: DETAILS OF THE QUARK DYSON SCHWINGER EQUATION 

The quark Dyson-Schwinger equation of the present 2-loop model, Eq. (jH)), is 

S-\p) = ii> + m + j dk A^,(p - k) |^7m'^(^)7. + ^l,S{k)^,S{k)^,S{p)^}^ . (Al) 

Projection of Eq. (jAlll onto the basis of Dirac vector and scalar matrices produces coupled 
non-linear equations for the propagator amplitudes A and B introduced in Eq. ()19|) . With 
Lorentz invariants x = p^, y = k"^, z = k-p/ \/Wp^, and integration measure J dk = 
1 / {2it)^ J^ dy y j\ dz ^/l — z^, the equation then takes the form, given in Eqs. (PUJI and 
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(EU, 



A{x) 
B(x) 



hix) + 

Mx) + 



A{x)l2{x) + B{x)h{x) 

xA{xy + B{x)^ 
A{x) J2{x) + B{x) J3(x) 



(A2) 
(A3) 



xA{xy + B{x)^ 

where the Ii{x), Ji{x) are scalar integrals over non-linear combinations of A{y) and B{y) for 
spacelike y. Because of the gaussian form Eq. (jH)) for the kernel ^^u{p — k) in the present 
model, the results of the angular integrations required for evaluation of /j(x) and Ji{x) 
can be ex pre ssed in closed form. To this end one needs the two confluent hypergeometric 
functions J3J|, 



Zl =/\rftexp{Mt-l)}(l-F 
Zl =j\dte^^{^i{t-l)}{l-e 



7riFi(l/2,l;-2^) 



TX 



Fi(3/2,3;-2/i) 



(A4) 
(A5) 



where /i = l^fxyjiJ^ . These hypergeometric functions can be numerically evaluated for gen- 
eral complex argument /i using their series expansion (|yu| < 1), direct numerical integration 
(1 < |/i| < 17), or asymptotic formulae (|/i| > 17) [34]. There are two important cases: real 
/i and purely imaginary /x, whereupon the angular integrals reduce to modified Bessel or 
Bessel functions, respectively. For real /i (i.e., real x = p^ > 0) 



Zl = 7re-^/o(/x) 



TT 



Zl = -e~'^lM 



/^ 



(A6) 



whereas for purely imaginary /x (i.e., real x = p"^ < 0) 



7re^^Jo(^') 



TT 






(A7) 



where fi = zfi'. The quantities Ii{x),Ji{x) can then be expressed as the one-dimensional 
integrals 



h(x) 



,^ 2D p / {^/x-^yT] , . 



TCU^ JO 
2D r^ 



LU^ 



y + 2uj'^ 



UJ 



+ zn-y 2 + - 



UJ 



1 + 2^^ 



h{x) = G—^ f 

-RijJ^ JO 



ayyexp :r^ I x 



LU^ 



o'M 



z^. 



2 , 4 , 2/^+4^ 

. yo; +u} + 

^^4 \ X 



Z'^ 



2 2.4 2y^' + 4o^' 
-xy — byuo — y — 2ul) 

X 
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+^l{y) 



u 



iX . 






hix) 



2D foo 

TTUJ^ Jo \ 



X 



vyy 



LU^ 



\(^v{y)(Ts{y) X 



X 



71I L ^ 2^y^' + ^^' 



X 



JAx) 



Mx) 



2D 



dyyexp 



m + 

TCtU^ Jo 



{V^-y/yf 



O". 



UJ^ 



■iy) X 



G 



2D /•°° 



TTUJ^ Jo 



dy y exp 



iv^-^y 



LU^ 



\(^v{y)(rs{y) X 



ii 



Zj- (-xcu^ - yu^ - 3u') + Zl- [2xy + xu^ + yu^ + 3u^) 



Mx] 



<yl{y) 



^2D foo / (v/i-Vy)^ \ 
G — - I dyyexp ^r- I x 



TTUJ^ Jo 



UJ^ 



+ ^s(l/) 



Zl^ + Zlli-x-y-2u^) 



(A8) 



The equations can be solved numerically by iteration subject to the boundary conditions 
A{x) -^ 1 and B{x) -^ m for large spacelike x. Alternatively, once the spacelike integrals 
for Ii{x),Ji{x) have been obtained, one may seek to utilize the polynomial structure of 
the equations in the explicitly appearing amplitudes A{x) and B{x). For completeness, we 
outline the derivation of the polynomial form. 

Elimination of either A or B from the right-hand side of Eqs. ()A2|1 and ()A3|1 produces 



xA^ + B^ [{A - h)J, -{B- Ji)h] - (/2J3 - J2h)A 



xA^ + B' 



[{B - Ji)/2 -{A- Ji) J2] - (/2J3 - ^2/3)5 








(A9) 
(AlO) 

(All) 



where we have dropped the argument x from A and B. This can also be written as 

{I2J3 - J2/3) _ iA- h)Js -jB- J{)h _ {B- Ji)h -{A- h)J2 
[xA2 + 52] A B 

The second equality can be multiplied out to give an expression for B^: 

B^h = B [A(J3 - h) + J1/3 - hh] + A^J^ + A{J,h - hJ2) . (A12) 

This expression is now used to eliminate any factors B^ in Eqs. ()A9|) and ()A10|) . yielding 

B [A{E -C)- hE] + A^D + A{JiC - hD) - hE = , 
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B 



A\J^E + h{C- 2E)) + AE{2hh - 2hJs + Ls-Ji) + ^(AVs - hh-Ji - /| 
+A\J2iE -C) + A\C{hJ2 - Jih) + E{JJ2 - 2/1J2)) + AEhihJ^ - Jih) = 



where C = x/| + /|, -D = X/3J3 + I2J2 and E = I2J3 — J2h- These last two equations 
give two equivalent expressions for B so by eliminating B one gets a polynomial equation 
involving only A: 

=A^ [Cxh + DJ3 - E{2xh + ^2)] 

+^3 [-Cx/1/3 + D{-3IiJ3 + 2/3J1) + E{4xlih + 3/1 J2)] 
+A2 [C(-/2/3 + I3J1) + D{-4hhJi + 3/1V3 - /|) 
+E(-2a;/i2/3 - 3/1V2 + 2/2/3 - /3^3 



+A 



-C{hhJl + ^'^i) + Dihll + 2/i'/3 Ji - /1V3) 
+E{-2hl2h + 2/1/3 J3 + /1V2 - Ipi] 
+E{hllJ^-llh.h + ll) , (A13) 

from whose solutions one can construct B using 

A^/^ + A{J,C - /i/}) - /3/^ 
^- AiE^cY^hE • ^^^^^ 

The fourth order polynomial Eq. ()A13|) has in general four solutions Ai, . . . , ^4 at each point 
X with an associated /?i, . . . , /?4. The boundary conditions of perturbative behavior in the 
asymptotic spacelike region (x -^ 00), together with the criteria of continuity specify the 
physical solution. 
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